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APPLICATIONS TO AERONAUTICS OF THE THEORY OF TRANSFORMATIONS OF NONLINEAR SYSTEMS 
George Meyer,* Renjeng Su,** and L. R. Hunt*^ 

NASA Ames Research Center, Moffett Field, California, U.S.A. 


ABSTRACT 

We discuss the development of a theory, its application to the control design of nonlinear sys- 
tems, and results concerning the use of this design technique for automatic flight control of air- 
craft. The theory examines the transformation of nonlinear systems to linear systems. We show how 
to apply this in practice, in particular, the tracking of linear models by nonlinear plants. Results 
of manned simulation are also presented. 


INTRODUCTION 

Suppose we model a physical plant by a nonlinear system 

m 

x(t) = f(x(t)] + ^ ^ u^(t)g^[x(t)] , (1) 

i=l 

where f, gj , ..., g^ are c“ vector fields on IR^ and f(0) = 0. If we are to have the output of 
this plant follow a particular path, then we have a difficult problem to consider. However, if there 
are new state space coordinates and new controls under which equation (1) becomes a linear system, 
then our task appears to be much easier because of the known results for controller design on linear 
systems . 

We feel that the following problems are thus of interest: 

(a) Find necessary and sufficient conditions for the system (1) to be transformable to 
a controllable linear systems. 

(b) Show how to use these transformations so that the controller design for nonlinear systems 
can be reduced to that of linear systems. 

(c) Apply the above theory to the field of aeronautics. 

In the next three sections of this paper we discuss tlie solutions of these problems. 


TRANSFORMATION THEORY 

The cl assi f ica tion of those nonlinear systems that can be transformed to linear systems is actu- 
ally a subproblem of a much deeper result, the construction of canonical forms for nonlinear systems. 
We are presently developing a theory for such canonical forms, and in the case that a nonlinear sys- 
tem is transformable as in this paper, the canonical form is actually the Brunovsky (ref. 1) form for 
a linear system. 

Here we concentrate on the transformation theory developed in references 2, 3 and 4. Other 
significant research in this area is due to Krener (ref. 5), Brockett (ref. 6), Jakubcyzk and 
Respondek (ref. 7), and Hermann (Tlie Theory of Equivalence of Pfaffian Systems and Input Systems 
under Feedback). We also refer to the early work of the first author in references 8 and 9. 

If we are to map our nonlinear system (1) to a controllable linear system, we may as 

well assume that this linear system is in Brunovksy (ref. 1) canonical form with Kronecker indices 

m 

iCj, < 2 * •••» satisfying 52*^^ ~ ^ ^rid ic^ > <2 - • • • - Hence this system is 

y = Ay + Bw , (2) 

where A is nxn,Bis nxm, w= (w^ jW^ , ..., w^) are the new controls, A is equal to 
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and B is given as 



The transformatio': results we present are actually local (in some open neighborhood of the 
origin in (x], X 2 , space), and global theorems are found in reference 3. We simplify nota- 

tion by saying jpr' when we actually mean an open neighborhood of (0, 0, ...» 0) in ]R^. However, 

F”' means (uj , U 2 , u^) space (or (wj, W 2 , w ) space) and this is not local. 
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We discuss the allowable transformations mapping system (1) to system (2). We want a c“ map 

Y “ (yi » 72 j •••» yn» ^^2’ •••» '^m) “upping R” x R™[(xj, X 2 » x^, uj , U 2 , .... u^) space] 

to IPP X [(yj, yj, .... y^, Wj, W 2 , ...» Wjjj) space] that satisfies the following conditions: 

J 

1. Y maps the origin to the origin, 

2* yi> 72» •••» 7n functions of xj , X 2 , . . . , Xj, only and have a nonsingular Jacobian 

matrix, 

3. wj, W 2 , . . . , Wjj^ are functions of xj , X 2 , ..., x^, uj, U 2 , ..., Uj^ and for fixed , X 2 , 
Xn, the m X m Jacobian matrix of , W 2 , ..., Wjjj with respect to ui , U 2 , . . - , u^j 
is nonsingular, 

A. Y maps system (1) to system (2), 

5. Y is a one-to-one map of ^ R onto R x R , 

Next we introduce some basic definitions from differential geometry. 

If f and g are c“ vector fields on R^ , the Lie bracket of f and g is 

•^*8] = If f - U 8 


where "If and are Jacobian matrices. We let 

(ad°f,g) = g 
(ad*f,g) = [f,g] 

(ad^f.g) = (f,(f,g)] 

(ad*^f,g) = lf(ad’''^f,g)] . 

A collection of C“ vector fields hj, I12 , h^ is involutive if there exists C~ func- 

tions 7 . such that 
1 J X. 

r 

(h.,hj](x) = y ] 1 i i> j - i j • 

k=l 


Lt't (•»•) denote the duality between one forms and vector fields. If u' = wj dX| + lO^dx^ + 
... + (ii„dx is a differentiable one form and f a vector field on R , then 

n n 

(w,f> = uif] + 0)2 f? + ••• + 

To state the main result from reference 4 giving necessary and sufficient conditions for trans- 
forming system (1) to system (2) we need the following sets: 


C = {gj. 

(f.gll.-. 

• .(ad‘'*”'f,gi) 

»f?2 * I f 1 


, ,(ad' ‘"*f,g2 


[f.gj... 

. ,(ad"^'"”’f ,gj^) 

1 



c, ={g:. 

[f.gil... 

.,(ad’.i “f.gi) 

.g:.[f, 

«g2 ] » • • • 

, ,(nd' ) 'f.gr. 


[f.gj... 

.,(ad’'j"'f,g|jj) 

1 for 

j=l,2,. 

, . . ,m 


Theorem 2. 1 There exists a transformation Y = (yj , y^y , ..., satisfying condi- 

tions i) through v) above if and only if on 


1 ) 

2 ) 

3) 

Let 
St rue ted 


tlie set C 
each set 
the span of 

oi = <iy 02 
in reference 


spans an n dimen.sional space, 

is involutive for j = l,2,...,m, and, 

Cj equals the span of Cj ^ C for j = l,2,...,m. 

= Kj + ^ = *^1 ^'2 + '*■ ^ Then the 

4 by solving the partial differential equations 


transformation 


is con- 


<dyj,(ad^f,gj)) = 0,j=0,l 
<dyo^+j.(adJf,gp> = 0,j=0,l 


,C|-2 and i=l,2,...,m, 
,k;2“2 and i = l,2,...,m. 
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( 3 ) 


^ dy^ , (ad^ f , g^)) ** 0* and 1^1,2) • • • ,m, 

1D“1 

TTl 

i=l 

m 

i=l 


m 
i-1 

where the matrix 

<dyi, (ad'^^"*f,gj)> . . 

^'^yoj+p (ad’^2-lf gj)) ^ ^ 

<dy ^.,, (ad*^™ *f,gi» . . 

m -1 

is nonsingular. 

It can be shown that matrix (A) being invertible means wc can solve for u^, , ...» Uj^ 

in terms of wj , w:> , . . . , Wj^ in the last m equations in equation (3). 

Equation (3) can be formally solved by considering a sequence of ordinary differential equations 
as in reference A, but we shall not mention details here. 

If a nonlinear system is transformable to a linear system, we study the process of using tin* 
transformation to construct a controller for the nonlinear system. 


<dyi, (ad'^*”*f,g^)> 

<dy^^^,.(ad‘'-2-4.gJ>] 


• <dy ,(ad'^"> 'f,g 


» 


(O 


TRANS F0RM,\T IONS IN CONTROLLER DESIGN 

Let Y = ^2* •••♦ ^’n» wj, W 2 , ...» w^) be the transformation from system (1) to system (2) 

as before. The structure of the control system using transformation theory is illustrated in 
figure 1. The design scheme is implemented on the "linear part" of the diagram, and this system 
is in Brunovsky form, 

U’e ask that the output of the nonlinear system follow a particular path which corresponds to a 
trajectory for the output of the linear model. If we know how to design for the linear system, tlion 
we actually have a tracking of a linear model by a nonlinear plant. 

Linear design is used to generate an open loop command w^ , for the system (2), and we find 
the corresponding y coordinates y^, by plugging w^. into equation (2). Tlic transformation Y 
maps the measured x space to y space and y is compared to and the difference is an 

error Cy. The regulator yields a control 6w which sends ey to zero, and variations in plant 
dynamics' and disturbances are compensated for in this way. 

The controls w^ and 6w are added and transformed through the inverse map R (actually 
w^ + 6w is substituted into the last m equations in equation (3) and u = (uj, U 2 , .•>, Uj^) is 
generated) to obtain a control which is applied to the plant. Thus we have an exact model follower, 

■ ::k 1 the difficult problem of finding an open loop control and the regulator control are constrained to 
tlir 1 incar system. 

The remainder of this paper contains the application of the transformation theory to aeronautics. 


AUTOMATIC FLIGHT CONTROLLER DESIGN 

The aircraft will be represented by a rigid body moving in 3-dimensional space in response to 
gravity, aerodynamics and propulsion. 
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The state 


X 



e XC X ]R3 X S0(3) X ]R3 


(5) 


where r and v are inertial coordinates of body center of mass position and velocity, 
C is the direction cosine matrix of body fixed axes relative to the runway fixed axes 
and 0 ) is the angular velocity. 


respectively; 
(inertial) , 


The controls, 


u 



U C X IR 


( 6 ) 


is the 3 -axls moment control such as ailerons eUyator and rudder in a conventional air- 


tlu or'roU Sclict’^^l “clira^rtl^rV^torroUe^i^ri^'a^helicopt and u^ controls 

power - throttle in a conventional aircraft, and the main rotor collective in ® 

state eq^tion consists of the translational and rotational kinematic and dynamic equations. 


V 


V * f^(x,u) 
c = S(ti))C 


(7) 




wherc^ and T* are the total force and moment generation processes and xcX. We wish to tran..form 
equation ( 7 ) into a linear system. 


In general, is invertible with respect 

cla.ss of helicopter maneuvers being considered 
to tho (scalar) pair Thus, a function 


to the (vector) pair (u,u‘*), and, for the specific 
(i.e., no 360° rolls), f*" l.s Invertible with respect 
h:X U can be constructed such that if 



h(r,v,C,u,u)OiV3o) 


( 8 ) 


then 


UQ 


'30 


(9) 


for all admissible maneuvers. That is, ans«lar and vertical accelerations can be ^hosen as the 
.set of independent controls in which case the state equation may be written as follows 


V » f°(r,v,C,V 3 o) + cf'(r,v,C,V 3 o,u,uio) 
C - s(u))C 


(in) 


0) = (jJO 


where r = 1, and f ' (r ,v,C, 0,0,0) = 0 for all admissible maneuvers. 

The function is invertible with respect to ((vi ,vz .E 3 ( 0 )) .C) where £ 3(0 is an elementary 

rotation about the z-axis and represents the heading of the helicopter. Thus, a function h :F 
X SO(2) * S0(3) can be constructed such that if 


Co » hf (r,v,vo,E3(ito)) 


( 11 ) 


then 


Vo 


( 12 ) 


constructed as follows 
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( 13 ) 


Co = h^(ro,vy,v(t),E3, (\i»(t)) 

(j 0 “ ^ (cc ) 

(Lq = 

uq = h(ro, VqjCq, 030.^0,0) 

where the function q extracts (d from CC^ = S(w) . The required time derivatives in equation (13) 
can be computed provided that the path (r,E 3 ) i| generated by the system diagrammed in figure 2 
where • represents- a scalar Integrator and y^ the control (w in the previous section is y^ here). 

We construct an approximation to the linearizing transformation as follows: Yi, Ri , Q are con- 

structed so that 


y = Y(x) = Yo(xo) + Yi 6 x = 

u « R(x,y^) * uq + Ri 6 y^ + 

Here Y^ ^is the transformation when e = 0, and 6 x (and 
xr. (and yp given in equation (13) (and figure 2). 

From equation (10) with c « 0, it follows that (C = 

( 6 r)* = 6 v 


yo + Yi5x 

Q5x 


(14) 


6 y^) is the perturbation about the nominal 


(I + S(e)Co) 


( 6 v)- 


iL 

3r 


5r 


0 



6 v 


0 



c 



(13) 


(t)’ = 6 u) 

(5u)) = 5 liq , 

vliero t is attitude perturbation. 

Thf pattern i>f equation (15) after some rearrangement of coordinates is sliown in equation (16). 


r;- 

V] 

V 2 

*•2 
•-3 
^3 
10 1 
(02 

^‘3 

V3 


0 

I 

— 

0 



0 

0 

^1 

C 2 

C 3 

0 


0 

0 

0 

I 

0 

0 

0 

0 


T 2 

vi 

V 2 

^-2 

^‘3 

^■1 

L >1 

(02 

(03 

/3, 


I 



(16) 


Tn the present case of the helicopter, Cj, € 4 , and C 5 are negligible. Their effect will be con- 
trolled by the regulator. 





The transformations 



( 17 ) 


(18) 


take the system in equation (16) (with Cj, C., , C*, = 0) into the canonic system: 


(19) 


Thus, the linearizing transformation (Y and R in figure 1) is constructed. 

That the accuracy of the transformation Is adequate may be seen from the results of the simula- 
tion of the flight experiment to be briefly summarized next. 

The test consists in automatically flying a trajectory which exercises the system over a wide 
range of flight conditions as shown in figures 4 and 5. 

Thus, the test takes the helicopter from hover (WPl) to high speed (150 ft/sec) turning accelera- 
tion, ascending flight. 

Figure 6 shows the resulting tracking errors. 

As can he seen, position tracking error e is quite small. The acceleration errors e^ , which 
is due to the neglected terms in the construction of the linearizing transformation is also quite 
small. In summary, the resulting performance of the system is good. 
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LINEAR SYSTEM REGULATOR TRANS- NONLINEAR 

FORMATION PLANT 


AND 

INVERSE 


Figure 1. Structure of the Control System 
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Figure 3. Experimental Flightpath Shown in Horizontal 
and Vertical Planes 
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